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Abstract 

F. Rohrlich has recently published two papers advocating a particular 
delay-differential (DD) equation as an approximate equation of motion 
for classical charged particles, which he characterizes as providing a "fully 
acceptable classical electrodynamics". We study the behavior in the re- 
mote past and future of solutions of this equation for the special case in 
which the motion is in one spatial dimension. 

We show that if an external force is applied for a finite time, some solu- 
tions exhibit the property of "preacceleration" , meaning that the particle 
accelerates before the force is applied, but that there do exist solutions 
without preacceleration. However, most solutions without preaccelera- 
tion exhibit "postacceleration" into the infinite future, meaning that the 
particle accelerates after the force is removed. Some may regard such 
behavior as sufficiently "unphysical" to rule out the equation. 

More encouragingly, analogs of the pathological "runaway" solutions of 
the Lorentz-Dirac equation do not occur for solutions of the DD equation. 
We show that when the external force eventually vanishes, the proper 
acceleration vanishes asymptotically in the future, and the coordinate 
velocity becomes asymptotically constant. 



1 Introduction 

The correct equation to describe the motion of a charged particle in flat space- 
time (Minkowski space) has long been a matter of speculation and controversy. 
The most-mentioned candidate has been the Lorentz-Dirac equation, written 
here for units in which light has unit velocity and metric tensor of signature 
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This is written in traditional tensor notation with repeated indices summed and 
emphasized in Greek; u — u' 1 denotes the particle's four-velocity, m and q its 
mass and charge, respectively, t its proper time, and Fjj the (antisymmetric) 
tensor describing an external electromagnetic field driving the motion. 

However, many objections have been raised to this equation. Among them 
are the existence of "runaway" solutions in which the acceleration increases 
exponentially with proper time even when the external field asymptotically 
vanishes. In some physically reasonable situations, all solutions are runaway 
(Eliczer, 1943; for later references see Parrott, 1987, Section 5.5). Even in fa- 
vorable cases in which non-runaway solutions exist, they may exhibit so-called 
"preacceleration" in which the particle begins to accelerate before the external 
field is applied (Rohrlich, 1965). 

F. Rohrlich (1997, 1999) has recently advocated a new equation of motion 
which (Rohrlich, 1999) claims without proof "[has] no pathological solutions". 
This new equation is a delay-differential equation given below, which we shall 
call the "DD equation" . The present paper derives from a study of the DD 
equation with the aim of determining the correctness of such claims. 

(Rohrlich, 1997) claims that the DD equation has neither preaccelerative nor 
"runaway" solutions. We shall show that for the case of a nonzero external force 
applied for only a finite time, the DD equation does admit preaccelerative solu- 
tions, and that these are all "runaway" in the weak sense that the acceleration 
does not vanish asymptotically in the distant past as one would expect. How- 
ever, this does not rule out the DD equation as a realistic equation of motion 
because we also show that these preaccelerative solutions can be eliminated by 
appropriate choice of generalized initial conditions (defined below). But then 
another problem arises: assuming this choice, most solutions exhibit "postaccel- 
eration" , meaning that the acceleration persists after the external field is turned 



Postacceleration is not as bad as preacceleration because there is no violation 
of causality, but postaccelerative solutions could be considered pathological. 
Suppose we are sitting in a room shielded from electromagnetic fields watching 
a beam of identical charged particles shoot in the window. It might seem strange 
if some of the particles speeded up, while others slowed down, for no apparent 
reason, according to their past histories. This is what the DD equation predicts. 

(Rohrlich, 1997) presents the DD equation as a modification of the following 
equation proposed by (Caldirola, 1956): 



This is written in units in which the velocity of light is c, and t := 4<j 2 / (3mc 3 ) w 
1.2 x 10~ 23 sec is a constant with the dimensions of time. With appropriate units 
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-[«'(T - T ) - U a (T - T )u a (T)u l (T)} + ^^(t) = 
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employing our convention that the velocity of light is unity, this can be written: 

-[u l (T-T Q )-U a (T-T )u a (T)u*(T)]+qF\u a =0 . (3) 
TO 

The motivation presented by (Caldirola, 1956) involves starting with the 
Lorentz equation (the classical force equation ignoring radiation reaction), 

m^—=qF\u a , (4) 
dr 

and attempting to replace du l /dr by the difference quotient 

u 1 {t) ~ u 1 (t - T ) 



TO 



(5) 



But since both sides of are orthogonal to u, while JSJ need not be, one ought 
to project (JSJ) into the subspace orthogonal to u{t). This projection will kill any 
multiple of u(r), so the resulting equation can be rewritten as 

Rohrlich (1997) suggested the following modification of J3J. This is what we 
are calling the DD equation: 

dn % 

m 1 — = f l (T)+m 2 [u l (T-T 1 )-u a (r-T 1 )u a (Ty(T)] . (6) 

Here /(r) = f l (r) is a four- force orthogonal to u(t), mj and m-i are presum- 
ably nonzero parameters associated with his motivation of the right side as an 
approximation to the self-force on a spherical surface charge, and n is a positive 
parameter. (The sign of the second term in brackets differs from his because his 
metric is opposite in sign to ours.) To avoid degenerate cases, we will assume 
below that m 2 and mi are nonzero; actually, they would be expected to be 
positive. 

(Rohrlich, 1997) attributes the conjecture of © to (Caldirola, 1956), but JHJl 
is not mathematically equivalent to ©, and (Caldirola, 1956) did not propose 
it. For example, Caldirola's equation for a compactly supported force admits 
periodic solutions (to which Caldirola attached great importance as possibly 
describing spin-like internal particle motions), while the DD equation does not 
(as the analysis of this paper will make clear). Note also that the DD equation 
P is a delay-differential equation, while Caldirola's equation is a difference 
equation involving no derivatives. 

The abstract of (Rohrlich, 1999), begins as follows: 

"The self-force for the classical dynamics of finite size particles is 
obtained. It is to replace the one of von Laue type obtained for 
point particles. ..." 1 

(Rohrlich, 1999) emphasizes the contribution of (Yaghjian, 1992, Appendix D., 
Equation D.19). Yaghjian derived a rest-frame version of the DD equation as 

lr The equation "of von Laue type" is the Lorentz Dirac equation tTl . 
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an uncontrolled approximation 2 to the self-force on a charged particle modeled 
as a "rigid" sphere. 3 

Yaghjian's equation is written for a particle slowly moving relative to some 
given Lorentz frame. If the given Lorentz frame is the particle's rest frame at a 
given proper time, and it is written in four-dimensional notation, then it becomes 
the DD equation. Thus the DD equation © is a relativistic generalization of 
Yaghjian's equation, and its statement in (Rohrlich, 1997) appears to be the first 
explicit statement in the literature. Moreover, (Rohrlich, 1997), goes beyond its 
mere statement to present it as a "fully acceptable classical electrodynamics". 
By comparison, Yaghjian obtained it as an approximation without comment on 
the approximation's domain of validity. 

Because of this and the above quote from the abstract of (Rohrlich, 1999) 
stating that it "obtained" the equation, an earlier version of the present paper 
called the DD equation © "Rohrlich's equation" . However, when Professor 
Rohrlich submitted a report as an identified referee for that version, he strongly 
objected to the name "Rohrlich's equation" , attributing it instead to Caldirola 
and Yaghjian. In deference to his wishes, we sidestep the convoluted question 
of the origin of equation JBJ by employing the neutral name "DD" equation. 

Although Rohrlich's motivation for the DD equation |JB} involves thinking of 
the particle as a charged sphere, the equation of motion itself is the equation for a 
moving point, the center of the sphere. Thus the DD equation is mathematically 
the equation of a point particle, and we shall consider it as such; the motivation 
of the charged sphere will not enter into our considerations. 

To study jnjl, it will be convenient to eliminate inessential constants by in- 
troducing a new time unit equal to t\ old time units and by using available 
notational freedom to absorb the two parameters mi and into a single pa- 
rameter m, obtaining 

du* 

m— = f + [u l (T-\)-u a {T-l)u a {r)u\T)] . (7) 

CLT 

This can be done by dividing both sides of by m^, and renaming / and m\. 

The resulting equation 0, which we shall call the "normalized DD equa- 
tion" , is equivalent to the original DD equation ||BJ| in the sense that knowledge 
of all solutions of either translates immediately into knowledge of all solutions of 
the other. In response to a referee's skepticism about this, we added Appendix 
2, which works out in detail the precise relations between the two equations. 
These will also be briefly sketched just below for the reader's convenience. 

The analysis of the paper could be written without introducing the normal- 
ized DD equation J7J), and the conclusions would be essentially the same, apart 

2 For example, the derivation ignores nonlinear terms. "Uncontrolled" describes the fact 
that no attempt is made to quantitatively assess the magnitude of possible errors which such 
simplifications might introduce.) 

3 Yahgjian mentions that essentially the same equation was "stated without proof" by 
(Page, 1918), and gives references to other precursors. Also, Caldirola refers to precursors 
of his proposed equation. The history of these ideas is extremely tangled, and we make no 
attempt to present all aspects of it here. 
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from trivial changes of language. We chose not to do that because retaining 
the original, more complicated notation would only distract attention from the 
more fundamental ideas which we want to emphasize. 

To avoid any misunderstanding, we make explicit that the unknown function 
"u" in (7J is technically not the same as the "it" in JBJ, but the two "it"s 
are related in a mathematically straightforward and physically transparent way 
given below in equation J3J), so that conclusions about solutions of (JTJl can be 
directly translated into conclusions about solutions of (JSJ, and vice versa. 

More explictly, changing to the new time unit corresponds to multiplying 
the old metric by . This results in a new proper time f :— t\t and a new 
four-velocity u satisfying 

u(f) = t iU {t) . (8) 
Equivalently, for any real number s, 

u{s) = tiu(txs) . (9) 

This simple relation allows us to translate any conclusion about u into a corre- 
sponding conclusion about it, and vice versa. 

The original DD equation © is thus equivalent to an equation like Q but 
with symbols renamed: 

m^ = r + [u\f-l)-u a {f-l)u a {f)u\f)} . (10) 

(XT 

The details of the relations of the renamed symbols to the originals are given 
in Appendix 2. 

In the rest of the paper, we omit the bars for typographical simplicity, and 
analyze (JJJ) instead of QlOjl. leaving to the reader the straightforward translation 
of mathematical results about J7J to conclusions about ||BJ|. All of our con- 
clusions about IHl will be essentially identical to the corresponding conclusions 
about 1|1(J [I . apart from trivial changes in language or notation. For example, 
Theorem 5's conclusion that for eventually vanishing force, \im f ^ du/d? = 0, 
translates into the conclusion that solutions it of © satisfy lim r ^ du/dr = 0. 



2 The DD equation for one-dimensional motion 

We shall study the (normalized) DD equation Q for the case of motion in one 
space dimension. (Henceforth we omit the reminder "normalized".) Analysis of 
this special case turns out to be sufficient to answer the questions motivating 
this work. 

For this case, we may work in a two-dimensional Minkowski space with 
typical vector x = (a; , a; 1 ) and metric tensor g(x,x) — (x ) 2 — (x 1 ) 2 - Since 
u a u a = g(u, it) = 1, we may write 

u = (cosh 8, sinhfl) , (11) 
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where this defines the "rapidity" parameter 0. Define 

w := (sinh0,cosh0) , (12) 

so that w is a spacelike unit vector orthogonal to u, and any vector orthogonal 
to u must be a multiple of w. 

In the DD equation (7J, all of du/dr, /, and the bracketed term are orthog- 
onal to u, and hence multiples of w. This can be seen explicitly for the left 
side: 

du dO 
m— — m—w 
dr dr 

It is natural to name 

the "scalar proper acceleration" . 
Write 

f = Ew , 

where this defines the scalar E. When f l = F l a u a is a Lorentz force, E is the 
electric field (nominally relative to the Lorentz frame corresponding to the basis 
u, w for two-dimensional Minkowski space, but actually relative to the Lorentz 
frame determined by any orthogonal basis, as is revealed by straightforward 
calculation) . 

The bracketed term of Q, being a multiple of w, is equal to its projection on 
w. The projection of an arbitrary vector v on w is —v a w a w, so the projection 
of the bracketed term of (JZJ) on w is w multiplied by the scalar 

-(cos!i(9(t - l)smh0(r) - sinh 0(t - 1) cosh 0(t) = sinh(0(r - 1) - 6>(r)) . 

Projecting the entire DD equation Q on w yields the equivalent scalar equation 

dO 

mA = m— = £ + sinh(0(T - 1) - 6{t)) . (14) 

We shall call this the variant of the DD equation for one-dimensional motion 
the "DDI equation". 

Equation Ijl4|l may be regarded as a delay-differential equations of the general 
form 

^=$(t,A(t) ) A(t-1)) , (15) 

with $ a given function. If we imagine E(t) as a given function of proper time 
r, then ifH|l is of the form (|To|l with $(r, r, s) :— E(t) + sinh(r - s).) 

In general, the situation is more complicated because is defined by (|11|) . 
with u 1 :— dz l /dr, where z 1 (t) represents the particle's wordline. Then E is 
usually not explicitly given as a function of proper time, but instead is given as a 
function E(z, u) — E(z(t), u(t)) of the Minkowski coordinates and four-velocity, 
with its proper time dependence acquired indirectly from the time-dependence 
of the latter. However, we may still regard the equation l|14|l as of the form l|15fl 
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by imagining solving 1141) for 9(t), which determines u(t), then (by integration) 
z(t), and finally E(z{t),u{t)). (The existence of the solution u{t) will be 
established in the next section.) Following common abuse of notation, we write 
E{t) in place of E(z(r), u(t)). 

This shows that 9 does satisfy some delay-differential equation of the form 
l|15|) . We shall show that this severely restricts the form of the function 9. For 
example, we'll show that if r i— > E(t) has compact support, then 9 must be 
bounded on any semi-infinite interval [to,oo). 

Similar remarks apply to the equation obtained by differentiating (|14|l : 

m^- = ^ + cosh(0(r - 1) - 9(t))[A(t - 1) - AM] . (16) 
dr dr 

This may also be regarded as of the form l|15|) if we imagine that we have 
already solved for 9. At first sight this may seem strange because if we have 
solved for 9, then we also have A := d9/dr, so there is no need to solve (|16f) 
for A. Nevertheless, the observation that A satisfies an equation of the form 
l|16(l is useful because it severely constrains A. For example, we'll show that it 
implies that for a force E(-) with compact support, linv^oo A(t) exists. Then 
combining this with the above-mentioned fact that 9 is bounded will imply that 
in fact liniT-^oo A(t) = 0. 



3 General delay-differential equations of form 
(USD 

This section reviews some simple and well-known facts about general delay- 
differential equation of the form (|15|) . When discussing such equations, we will 
always call them "delay-differential" equations, which will never be abbreviated. 
We reserve the term "DD equation" for the particular special case {7J) or one of 
its variants such as ©■ 

Let the function $ in the delay-differential equation i|15|) be C 1 (i.e., con- 
tinuously differentiable) . Suppose that r t— > A(r) satisfies this equation. If we 
regard A(r) as given on some interval [n — 1, n], then (|15f) becomes an ordinary 
differential equation for A on [n, n + 1] of the form 

^- = #(t,A(t)), 7i<r<n+l , (17) 
dr 

which is covered by the standard existence and unigueness theorems. 

This observation reveals the general structure of solutions of QBt - Choose 
A(-) to be an arbitrarily chosen C 1 function on any closed interval of length 1, 
say the interval [0, 1], subject to the consistency condition 

^(1) = $(1,A(1),A(0)), (18) 

where the derivative in l|18fl is understood as a derivative from the left. Such a 
specification of A on an interval of length 1 will be called a generalized initial 
condition. 
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Then (|15fl determines a unique solution A(r) on some maximal interval 1 < 
t < 5 with S < 2. We shall show below that for the equations of interest to us, 
namely (|rl|) and l(Po)l . we actually have 5 = 2, and A(-) satisfies the equation 
on [1,2], where A'(2) is understood as a derivative from the left. Iterating the 
process yields a solution A on [0, oo) whose values on [0, 1] can be an arbitrarily 
specified C 1 function satisfying the consistency condition l(T%|) . 

Iterating to the left to obtain a unique solution on (— oo, 0] determined by 
A restricted to [0, 1] involves inverting s i— » $(r + 1, A(r + 1), s) for fixed r. For 
equations lfT^|> or ijPol) . this is trivial. For example for (|H|) . given 9(t) defined 
for < t < 1, simply define 9 on the "preceding" interval — 1 < r < 

d8 

9(t) :=sinh _1 (m— (t + I) - E(t + 1)) + 6(t + 1), -1 < r < 0, (19) 
dr 

where the derivative d6 / dr(r + 1) is understood as a derivative from the right 
when r + 1 = 0. 



4 Special cases of the general delay-differential 
equation ( 1151) 

4.1 Mathematical preliminaries 

The delay-differential equation l|15|l relates the solution A(-) on an interval [a, a+ 
1] to the solution on the "preceding" interval [a— 1, a] . The following proposition 
shows that for a class of equations which includes l|14f> and (|16|l (the latter with 
9 regarded as given, a priori), the maximum of A on an interval [a, a + 1] is 
no greater than the maximum on the preceding interval [a — l,a]. Similarly 
the minimum of A on [a, a + 1] is no less than the minimum on the preceding 
interval. 

For an arbitrary C 1 real-valued function A on the real line, and arbitrary 
a < 0, define: 

M [Qi/3] (A) := max{A(r) a < r < (i } (20) 
m M (A) := min{A(r) | a < r < (3 } (21) 



Proposition 1 Let r, s \— * Q(t, s), — oo < r, s, < oo, be a C 1 function such that 

for each r, s h J7(t, s) is a strictly increasing function satisfying 0(r, 0) = 0. 
Let X be a solution of a delay- differential equation of the special form 

^ = n(r,A(r-l)-A(r)) . (22) 

Then for all a, 

M [a , a +i] < M [a _ ha] (23) 

™[a,a+l] > ™[a-l,a] • (24) 
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Proof: For notational simplicity we take a := 0. Thus we will prove that 

M [h2 ] < M M . 

The proof of the corresponding assertion for m, which only requires reversing a 
few inequalities, will be omitted. 

Let r max denote a point in [1, 2] with 

A(r roax ) = M[ 1)2 ] (A) . 

If Tmax is an interior point of [1,2], then 

= A (r max ) = n(r max , A(r max — 1) — A(r max )) , 

and since s t— > f2(r, s) is strictly increasing and zero only at s = 0, we have 

-^[o,i] > A(r max - 1) = A(r max ) = M[ 1]2 ] ■ 

If T max = 1, then 

M [0 ,i] > A(l) = A(r max ) = M [1>2] . 
If Tmax = 2, then we must have 

A' (2) > , (25) 

otherwise there would be points r < 2 with A(r) > A(2) (i.e., the graph of 
A would be going down at 2), contradicting A(2) = A(r max ) > A(r) for all 
1 < r < 2. Now from $B and E3, 

0< A'(2) = n(2,A(l)-A(2)) , 

so again using the fact that 17(2, s) is increasing with £1(2, 0) = 0, we conclude 
that A(l) - A(2) > 0. Finally, 

M [0;1] (A) > A(l) > A(2) = A(r max ) = M [li2] (A). 

I 



Now we prove that given an arbitrary C 1 specification of A(-) on the interval 
[0, 1] satisfying the consistency condition (|18f) . there exists a solution A to IL'l't 
defined on [0, oo) and taking the specified values on [0, 1]. By a C 1 function on 
[0, 1], we mean a continuously differentiable function, with the derivatives at 
and 1 understood as derivatives from the right and left, respectively. 

Proposition 2 Let t i— » ip(r) be a C 1 function on [0, 1] satisfying the consis- 
tency condition lltfy Consider the equation ]2%fy of Proposition 1: 

^ = 0(r,A(r-l)-A(r)) , (26) 

where satisfies the hypotheses of that Proposition. Then there exists a unique 
C 1 solution A defined on [0, oo) and satisfying A(r) = ip(r) for < r < 1. 
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Proof: If we regard A(r — 1) = tp(r — 1) as given on 1 < r < 2, then equation 
<|26|) is covered by the standard existence and uniqueness theorems for ordinary 
differential equations (e.g., Perko, 1996, Chapter 2). By these results, given the 
initial condition A(l) = ^(1), there exists a maximal interval [1,8), 1 < S < 2, 
such that there exists a C 1 solution A on [1, S) satisfying that initial condition. 
Moreover, if A(r) remains in a compact set for 1 < r < S, then 8 — 2 (Perko, 
1996, Section 2.4, Theorem 2). 

In other words, the only way the solution can fail to be globally defined is if 
it blows up. By Proposition^ our solution does not blow up; hence it is defined 
for < r < 2. Iteration produces a solution defined on [0, 00). I 

We can also iterate to the left to obtain a solution defined on (—00,00), 
assuming that for fixed r we can invert s 1— » f2(r, s). The inversion is trivial for 
equations of the form (|14(l and Ijl6|l . so we have: 

Proposition 3 Let ip be a C 1 function on [0, 1] satisfying the consistency con- 
dition UflJ) . Then for any C 1 function r 1— > E(t), there exists a C 1 solution 
t 1 ► 9(t) to the DDI equation \14\j , defined for —00 < t < 00 and satisfying 
9{t) = %1){t) for < r < 1. 

The same holds with A in place of 9 for the acceleration equation (|16[) , provided 
we regard 9 as given a priori and assume that E(-) is C 2 , so that dE/dr is C 1 . 

4.2 Preacceleration and solutions runaway in the past 

The simple propositions of the previous subsection imply quite a lot. Consider 
the DDI equation lfT4*|) 

d9 

mA = m— = E + sinh(9(T - 1) - 9(t)) . 
dr 

for a continuous force E applied for only a finite time, say for < r < Tf (i.e., 
E(t) vanishes off this interval). 

Suppose that A does not vanish identically on [—1,0]. Then A cannot vanish 
asymptotically in the past because by Proposition ^ 

M[_ li0] (A) < M[_ 2 ,-i] (A) < Af^s.-a] (A) < ... , 

and 

m[_ 1)0] (A) > m[_2,_i]04) > m[_ 3 - 2 \{A) > ... . 

The first equation tells us that if M[_ lj0 ](^4) > 0, then A assumes values at 
least as large as this infinitely often in the distant past. If, on the other hand, 
M[-i.q](A) < , then also m[_i^](A) < M[_ li0 ](A) < 0, so \A\ assumes values 
at least at least as large as m[_ li0 ](^4) infinitely often in the distant past. So 
the only case in which A could vanish asymptotically in the distant past is if 
m [-i.o](^4) — A/[_i.o](^4) = 0, but then A vanishes identically on [—1,0]. 
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Thus the only physically reasonable initial specification of A on [—1, 0] is 

A(t) =0 for -1 < t < . (27) 

A glance at <|l(i|) shows that this implies that A(t) vanishes identically for t < 0, 
and so is constant (i.e., the velocity is constant) for t < 0. 

So, we see that even for a compactly supported force (indeed, even for iden- 
tically zero force), there are solutions of the DD equation which do not vanish 
asymptotically in the past. These solutions also exhibit "preacceleration" . How- 
ever, for a compactly supported force, we can choose the initial specification l|27[l 
of A so as to eliminate these pathological solutions. 

The situation for a force which is not compactly supported seems unclear. It 
does not seem obvious how to choose the initial specification so as to force A to 
vanish asymptotically as r — > — oo when E does. Since the theory is physically 
incomplete unless one gives a prescription for an initial specification of 6 on 
some interval of length 1, this is a point which advocates of the DD equation 
should address. 

Suppose we have agreed on such a prescription. We might define the A de- 
termined by the prescription as preaccelerative if there exist two force functions 
Ei and E 2 with E 1 (t) = E 2 (t) for r < 0, but Ai(r) ^ A 2 {t) for some r < 0, 
where Ai and A 2 denote the corresponding accelerations obtained by solving 
the DD equation. It does not seem obvious that there is a prescription which 
will outlaw preaccelerative solutions. 

4.3 Postacceleration 

This section considers a force applied for only a finite time Tf, say < r < Tf. 
For such a force, a preaccelerative solution of the DD equation J7J or DDI equa- 
tion l|14l) will be defined as one such that the acceleration A of the corresponding 
solution does not vanish identically on (— oo, 0], and a postaccelerative solution 
is one whose acceleration does not vanish identically on [t/, oo). 

The discussion to follow explains why either preacceleration or "postaccel- 
eration" is essentially built into the DDI equation: one or the other must occur 
except in unusual special cases. This can be guessed from the observation that 
we have two ways to solve (uniquely) the DDI equation (Q): 

dO 

mA = m— -_ = E + sinh(0(r - 1) - 0(r)) . (28) 

We may solve (|28|l forward in time from a generalized initial condition u(t) = 
for — 1 < r < 0. We can also solve backward in time from a "generalized final 
condition" u(t) = for Tf < r < r/ + 1. There is no reason to imagine that 
these two solutions will coincide, and we shall now show that they almost never 
do. 

We may take tj — n, with n a positive integer. Consider the DDI equation 
l|28|) for Tf = n < r < n + 1. Suppose there is no postacceleration, so that 
A = dO/dr vanishes on this interval. Then 9 is constant on this interval, say 
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8(t) = k, a constant, for Tf = n < r < n + 1. Substituting these facts into l|28[) 
gives 

= + sinh(0(r - 1) - k) for n<T<n + l , 
which may be restated as: 

8(t) = k for n — 1 < r < n . 

Next consider equation (|28|l on the interval n — 1 < r < n. Using the above 
observation that 8(t) = k on this interval, which implies that A = on this 
interval, we see that 

E(t) = - sinh(0(r - 1) - k) for n - 1 < r < n . (29) 

This says that, given the values of 8(t) for n — 2 < t < n — 1 (i.e., the values 
of 6{t — 1) for n — 1 < r < n), there is a unique force function -E(t) on 
n — 1 < t < n which will eliminate postacceleration. Put differently, for any 
other force function on n — 1 < r < n, postacceleration mwsi occur. 

To complete the argument, suppose there is no preacceleration, i.e., A(t) = 
for t < 0. Then for some constant 8$, 0(r) — 8q for t < 0. Solve the DDI 
equation 114|) forward in time starting with this generalized initial condition, 
8(t) = 8q for —1 < r < 0. This will produce a unique solution 8(t) on the 
interval < r < n — 1, and this solution depends only on the values of E{t) 
for < r < n — 1. Substituting this solution in equation l|29(l shows that 
postacceleration must occur for any force function E(t) onn - I < t < n which 
does not satisfy equation l|29|) . i.e., for almost all force functions. 

The above argument becomes much simpler for the case n = 1, correspond- 
ing to a force applied only for a time interval of length 1, i.e., a time interval 
the length of the delay parameter in the DD equation. In Rohrlich's motivation 
of the DD equation, this time interval is the time light required to traverse the 
spherical particle's diameter, so it is probably too small to have any observa- 
tional interest, but an examination of this special case is a quick way for the 
reader to become convinced that it is mathematically unreasonable to expect to 
simultaneously eliminate both preacceleration and postacceleration in general. 

The simpler argument for n = 1 goes as follows. Observe as before that 
8(t) must be constant for n — 1 < 8 < n, i.e., 8{t) = k for < 9 < 1, and in 
particular, 8(0) = k. If there is no preacceleration, then 8(t) must be constant 
for t < 0, say 8(t) = 8q, for r < 0. In particular, 8q — 8(0) = k. Hence on 
[0,1], 

8( T - 1) = O = k for < t < 1 . 
Substitute this in Q29JI to obtain 

E(t) = - sinh(/c - k) = for < r < 1 . 

This says that for an force applied no longer than the time delay, either prec- 
celeration or postacceleration must occur except in the trivial case of identically 
vanishing force. 
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5 Forward asymptotics for eventually vanishing 
force 



This section proves that solutions of the DDI equation (|14|) for eventually van- 
ishing force have proper acceleration asymptotic to zero and asymptotically 
constant velocity in the future: 

lim A(t) = ; lim 6(t) = constant . (30) 

T — >OC T — ► OO 

Thus replacing the Lorentz-Dirac equation by the DD equation does eliminate 
the undesirable "runaway" solutions of the former for the special case of motion 
in one space dimension. Whether this is also true for motion in three space 
dimensions is unknown. 

Suppose that E eventually vanishes, and choose the origin of time so that 
E(t) = for r > 0. Then for < r < 1, A satisfies 

A A 

m— = cosh(6»(r - 1) - 9{t))(A(t - 1) - AM) (31) 
ar 

We shall consider the related equation 

^ = (f>M(gM - hM), < r < 1, with initial condition h(0) = g(l). (32) 
ar 

Here 4> will be a given continuous function on [0, 1], and g another function with 
the same domain. We will think of cf> as fixed until further notice, and of l|32|l 
as defining a mapping which assigns to each g another function h, namely the 
unique solution of l|32|l . We will show that this mapping is a strict contraction 
relative to a certain Banach space norm. 

Let C[0, 1] denote the real Banach space of all continuous real-valued func- 
tions on [0, 1] with the supremum norm || • Hoc,: for g £ C[0, 1], 

Halloo :=max{| 5 (r)| | < r < 1} . 

For g € C[0, 1], let M(g) := M[ 0i ii(g) and m(g) :— m[ 0j i]((?) denote the maxi- 
mum and minimum of g. 

Let B denote the Banach space which is the quotient of C[0, 1] by the one- 
dimensional subspace of constant functions. For g 6 C[0, 1], let g denote its 
image in B under the quotient map. Denoting by 1 the function constantly 
equal to 1, and by R the real field, the quotient norm || • || on B is: 

\\g\\i=M{\\g + al\\ 00 \a&Il} = ^[M{g)-m{ g )] . (33) 

The first equality is the definition of the quotient norm, and the last is a simple 
exercise. 

Let be a given function, and consider the linear mapping : C[0, 1] — > 
C[0, 1] defined as follows. For any g G C[0, 1], 

Q<P9 ■= h , (34) 
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where h is the unique solution of 

Define Q$ : B — * B to be the analog of on the quotient space B: for all 
ff6C[0,l], _ _ 

Q<p(g) : = Q<t>9 (35) 

We see that is well-defined by noting that for any a 6 R, 

d(/t j" r al) = § = 0(t)[<?(t) - ft(r)] = ^(r)[( 5 (r) + a ) - (ft(r) + a)} , 

and also (h + al)(0) = /i(0) + a = + a = (5 + al)(l). Thus if we alter g by 
adding a constant a to it, then we also alter the solution of (|32|l by the same 
additive constant, so that the h in l|3"5|) depends only on the equivalence class 
of g in C[0, 1] modulo the constants. 

Next we show that for any non-negative <f>, the mapping : B — > B is 
a strict contraction. Actually, we'll need the following stronger fact giving a 
uniform bound on \\Qc/,\\ for any bounded set of non- negative <fi: 

Lemma 4 For each non-negative <fi £ C[0, 1], \\Q</>\\ < 1. Moreover, given any 
constant ko, there is a constant k < 1 such that for all non-negative functions 
4> with II0H < kg, and for all g G C[0, 1], 

IIQ^II <k\\g\\ . (36) 
Proof: For each non-negative function <j> S C[0, 1], define 

ip(x) := / 0(s) ds . 
Jo 

Note that since <fi is non-negative, ^ is non-decreasing. 

Then for any g e C[0, 1], (Q^p)(r) = ft(r) with /i the solution to 

ft(r) = e^ (r )p(l) + f e^ s) (p(s)g(s) ds , 

Jo 

so 

h(r) < e-^g(l) + M(g)e-^) f e Hs)^L ds 



ds 

= e-^ T) g(l) + M(g){l-e-^ T '>) 
= M(g) + (g(l)-M(g))e-^ . 

Since 5(1) — M(g) < and ip is non-decreasing, 

M(h) < M(g) + [g(l)-M(g)]e-^ 

= M(g)(l-e^) + g(l)e-^ . (37) 

Similarly, 

h(r)>e-^g(l)+m(g)(l-e-^) , 
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and since <?(1) — m(g) > and ip is non-decreasing, 

m(h) > m{g) + [g(l) - mO)]e _ ^ (1) 

= m(g){l - e-^) + g(l)e-^ . (38) 

Subtracting the inequalities (|37J) and (|3H|> shows that for each non-negative <fi, 
Qrj, is a strict contraction: 

\\QMW = [M(h)-m(h)}/2 

< (l-e-^)[M(g)-m(g)]/2 = (l-e-^)\\~g\\. (39) 

Finally, the uniformity condition follows from l|39|) with k := 1 — e~ fc °, since 
11011 < k Q implies ^(1) < k Q and hence 1 - e - ^ 1 ) < 1 - e- fco . I 

Now we prove the main result of this section: 

Theorem 5 Suppose that the force r i— > E(t) in the DDI equation \14\j even- 
tually vanishes, meaning that there exists a proper time t$ such that E{t) = 
for t > To . Then there exists a constant 9^ such that 

lim A(t) = and lim 0(r) = ^ . (40) 

r — >OG T — 'OC 

Proof: We choose the origin of proper time so that E(t) = for r > 0. In 
solving the acceleration equation (|16|l , we assume that the solution A is already 
obtained for r < 0. Then the solution on [0, 1] is obtained by applying an 
opertor to the restriction A\[— 1, 0] of A to [—1, 0], translated one unit right. 

More precisely, define an operator T on C[0, 1] to be right translation by one 
unit: (T/)(r) := /(r - 1) for / e C[0, 1] . Let <f>{r) := cosh(6>(T - 1) - 6»(r)), 
< t < 1. Then 

for0<r<l, A(r) = Q^T(A|[-1,0])(t) . 

Though we don't know 0(r) for < r < 1, we do know from Proposition 1 
that it is bounded above and below by bounds no worse than on the preceding 
interval — 1 < r < 0. Hence we have the a priori bound 

Moo < fco := cosh(|M(0|[-l,O])| + |m(0|[-l,O])|) . 

The solution A(t) for 1 < t < 2 is similarly obtained, except that the Q$ is 
different because the <fr is different. However, Proposition ^ shows that we have 
the same a priori bound on the new ||^||oo< After n applications of Lemma 
we find that 

(Af(i4|[n-l,n])-m(A|[n-l,n]))/2= \\A\[n - l,n}\\ < k n \\A\[-l, 0]|| . 

This implies that A\[n — 1,".], becomes asymptotically constant as n becomes 
large. The constant must be zero because A — dO/dr, and the previous section 
showed that the rapidity 9 is bounded on [0, 00]. Hence 

lim A(t) = . (41) 
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Finally, we show that the rapidity (hence the velocity) becomes asymptoti- 
cally constant in the future: for some constant 9^, 

lim 0{t) = 6^ . (42) 

r — >oo 

From the last section, Mr nn+ -n(6) is a nonincreasing function of the positive 
integer n, and ni[ n n+1 ](9) is nondecreasing. This implies that if Q42[l were false, 
there would exist a positive constant i5 with M[ nn+ i^(9) — rri[ n n+1 }(9) > S for 
all n > 0. Choose points r+ and t~ in [n, n + 1] with #(r+) = M[ n n+1 ](8) and 
9{t~) — ni[ n ^ n+1 ](8). Suppose t~ < r+. Applying the Mean Value Theorem to 
the interval [t~,t+] yields a point r n between r~ and r+ with 

A t r S - g ( T ») - d ( T n ) > A Wl](g) - ™[n,n+l](g) g ^ 

Similarly, if r+ < r,7, then ^4(r„) < —6 for some t„ £ [n,n + 1}; both this and 
l|43|l contradict the previously established fact that lim T ^oo A(t) = 0.1 



6 Summary and assessment 

The above proofs were given in the context of the normalized DD equation 
J2J for motion in one space dimension and its equivalent formulation, the DDI 
equation l|14|) . Translating them into information about the original DD equa- 
tion JHJ discussed in (Rohrlich, 1997, 1999) yields the following observations, in 
which the motion is assumed restricted to one spatial dimension unless otherwise 
specified. 

1. An appropriate condition to guarantee a unique solution is a specification 
of the rapidity 9{t) on some closed interval [a, a + ri], subject to the con- 
sistency condition (|18|l (written there for a :— and t\ = 1). We refer to 
this specification as a "generalized initial condition" . 

There exists a unique C 1 solution 9(-) of (|14|) satisfying any C 1 generalized 
initial condition, and hence a unique four- velocity u(t) = (cosh 6{t), sinh 9(t)) 
and worldline r i— > z(r) with dz/dr = u(t) (by integration, with arbitrary 
specification of z(0 j). 

2. For the case of a force r E(t) applied for only a finite time (i.e., 
compactly supported force), there exist solutions which do not vanish 
asymptotically in the distant past. Such solutions seem unphysical. For 
simplicity of language, we will call them "past runaway" solutions even 
though we do not know that the acceleration becomes unbounded in the 
infinite past. 

However, for compactly supported force, there exists a unique choice of 
generalized initial condition which eliminates these past runaway solu- 
tions. This seems the physically relevant choice of initial condition. 
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3. For a force not compactly supported, it has not been proved that there 
is an appropriate choice of generalized initial condition to eliminate past 
runaways. Indeed, it seems unclear what such a choice might be. If the 
DD equation is to be used to predict the motion of actual particles, some 
definite procedure for specifying the generalized initial condition would be 
necessary. 

4. For a force E which eventually vanishes (i.e., for some To, E(t) = for 
t > To), the acceleration does converge to zero in the future, and the 
velocity becomes asymptotically constant. Thus for the DD equation for 
motion in one space dimension, analogs of the pathological "future run- 
away" solutions of the Lorentz-Dirac equation do not occur. It is unknown 
whether future runaways are also impossible for general three-dimensional 
motion obeying the DD equation. 

5. For compactly supported force, there exist "preaccelerative" solutions (de- 
fined as solutions with nonzero acceleration before the force is applied). 
Such solutions seem unphysical because they violate "causality". How- 
ever, the same generalized initial condition which rules out past runaways 
also rules out these undesirable preaccelerative solutions. 

6. For compactly supported nonzero force and generalized initial condition 
chosen to rule out preacceleration, the resulting solution necessarily ex- 
hibits the unusual property which we call "postacceleration" : the accelera- 
tion persists after the force is removed, and in fact into the infinite future. 
Some may consider this strange behavior sufficiently unphysical to rule 
out the DD equation. Others may welcome the unusual prediction as a 
potential physical test which if verified, would constitute striking evidence 
for the equation. 

Given that known motivations of the DD equation arc not fundamental (e.g., 
employ uncontrolled approximations obtained by ignoring nonlinear terms), and 
given the uncertainties noted above about the mathematical appropriateness of 
the equation, it seems that more work would be needed to convert it from a 
speculative proposal to an accepted physical principle. 

It is often said that a resolution of the logical problems of classical electro- 
dynamics must come from more fundamental physical principles of quantum 
mechanics. This may indeed turn out to be true, but it is not the only pos- 
sibility. It would seem strange if a logically consistent and physically sensible 
classical electrodynamics were inherently impossible. Although nearly all of 
physics rests in principle on quantum mechanics, there do exist consistent and 
sensible classical theories in many areas. The present theory of classical charged 
particles cannot be considered sensible. 

It is interesting that the Lorentz-Dirac equation has survived as the principal 
candidate for a classical equation of motion despite predictions so bizarre (e.g., 
Eliczer, 1943, Parrott, 2002) that no one will admit to believing them. The 
reason for the survival may be the fundamental nature of Dirac's derivation of 
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the equation. If one accepts mass renormalization (admittedly controversial), 
then one can convincingly obtain the Lorentz-Dirac equation from the principle 
of conservation of energy-momentum with no approximations whatever. One 
does not lightly discard such mathematically tight arguments. 

The motivations of many other proposed equations of motion are aesthet- 
ically less pleasing. Most consist of more or less ad hoc modifications of the 
Lorentz-Dirac equation which do not obviously lead to its bizarre predictions. 
(This is not the same as saying that the modifications obviously do not lead to 
similarly bizarre predictions!) 

Assuming that one is willing to believe in postacceleration, the DD equation 
seems physically possible, but not physically compelling. It is not clear in what 
sense, if any, solutions of the DD equation conserve energy-momentum. It would 
be desirable to find a convincing, mathematically rigorous way to relate the DD 
equation to the principle of conservation of energy-momentum. 

7 Appendix 1 

When the body of this paper was written, I was unaware of seminal work of 
Ryabov (Ryabov, 1960 , 1961, 1963, 1965) on delay differential equations, sub- 
sequently refined and extended by (Driver, 1968). This work does not directly 
apply to the DD equation because this equation does not satisfy Driver's hy- 
potheses. However, it has the same structure as equations considered by Driver. 

If sinh were a bounded function, then Driver's results would apply to Rohrlich's 
equation (|14fl . Thus we can hope that conclusions similar to Driver's might ap- 
ply to the DD equation. If so, they would shed light on some of the problems 
mentioned in the Summary and Assessment section. 

This appendix outlines Driver's conclusions in the context of the DD equa- 
tion. It is mathematically based on the excellent semi-expository paper (Driver, 
1968), but the presentation follows the introductory sections of (Chicone, 2003), 
which is based on Driver's work. The latter presents a point of view particularly 
congenial to a discussion of the DD equation. 

We shall be concerned with delay-differential equations of the form (15), but 
for easy exposition we need to reinsert the positive delay parameter t% (which 
was set to 1 in (15) by appropriate choice of units): 

^ = *(r,A(r),A(r-r 1 )) . (44) 

In equation (15), the function A was a real- valued function, but the case in 
which A takes values in R n is no more difficult to discuss, and the results that 
we shall present apply also to this case. Therefore, we assume that A(r) £ R n . 

The specialization of Driver's work of interest to us applies to equations of 
this form with $ continuous and bounded. The hypothesis that $ be continuous 
and bounded makes the results easy to state. It can be relaxed, but not so far 
as to admit the Rohrlich equation. More precisely, the hypothesis that $ be 
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bounded can be relaxed to an assumption that satisfy a certain Lipschitz 
condition. 

A very important hypothesis for the results of Ryabov and Driver to be 
described is that T\ be sufficiently small. That is why we needed to reinsert the 
delay as an explicit parameter in (|44H . 

Previously we normalized n to 1 by appropriate choice of units, and the 
reader may be puzzled as to why we can't maintain the same normalization 
here. We could, but it's not convenient. The reason is that normalizing t\ to 
1 correspondingly changes the Lipschitz constant, so that the restriction that 
T\ be sufficiently small is replaced in general by a restriction on the Lipshitz 
constant, or by a restriction on the bound of 4> under our simpler hypothesis 
that $ be continuous and bounded. From a physical point of view, it is more 
natural to phrase the hypothesis for the Ryabov/Driver results in terms of a 
restriction on the time delay instead of a restriction on the Lipshitz condition 
or the bound of $. 

A priori, a mathematically appropriate initial condition to specify a unique 
solution for (|44|l is a specification of the function A(-) on an interval [0, ti) 
subject to the consistency condition dX/dt{r\) = $(0, A(n), A(0)). Thus the 
solution space of H44fl is infinite dimensional. 

For the DDI equation (|14|l (corresponding to n = 1 in (|44l0 . the physical 
expectation is that there should be precisely one physically relevant solution 
for each ordinary initial condition A(0) = Ao € R 1 . Thus the space of mathe- 
matical solutions of (|44|l is too large — one needs to find a condition which will 
reduce the infinite-dimensional solution space to a one-dimensional space. The 
Ryabov/Driver theory suggests such a condition, as we shall now discuss. 

A solution A of (|44|l is called a special solution if 

sup e-l T l/ Tl ||A(r)|| < oo . (45) 

— 00<T<00 

Thus a special solution is one which does not increase at too fast an exponential 
rate. In particular, all solutions which are not special are "runaway" (either 
in the past or the future), and therefore presumably unphysical for physically 
possible forces. 

Driver's results imply that for every Ao £ R n , there exists a unique special 
solution (to with $ continuous and bounded) satisfying A(0) = Ao- Thus 
the space of special solutions is a manifold of dimension n, coordinatized by the 
map A i— > A(0). 

In addition, Driver showed that this manifold is an attractor in the sense that 
every solution approaches (is "attracted" to) a special solution (exponentially 
fast) in the infinite future. Thus the manifold of special solutions models the 
long-term behavior of all solutions. 

The DD equation (J7J in m-dimensional space ((m + l)-dimensional space- 
time), can be reformulated as an equivalent equation of form (|44|) with n = m. 
Here we are primarily interested in the case m = 3, and secondarily in m = 1. 
The reason that @ with m — 3 is not immediately of the form 144(1 is that a 
four- velocity u must satisfy the auxiliary condition u a u a = 1. Taking this into 
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account reduces it to an equivalent equation of form (|44l) with m — 3; we omit 
the details. Although the corresponding $ on its right side is not bounded, we 
can hope that the conclusions of the Ryabov/Driver theory might still apply. 

Suppose that the force in the DD equation Q is continuous and compactly 
supported (hence bounded), and consider the equivalent equation of form l|44|l 
with to = 3. Then all physically relevant solutions should be bounded, hence 
special (because all non-special solutions are "runaway", as noted above). If the 
conclusions of the Ryabov-Driver theory extend to this more general situation, 
then the space of physically relevant solutions has dimension no greater than 

TO. 

In general, special solutions can increase exponentially fast, so we cannot 
immediately identify special solutions with "physical" solutions. However, it is 
unknown whether exponential growth of special solutions can occur for the DD 
equation in four-dimensional spacetime. 

If there were exponentially increasing special solutions to the DD equation 
in four-dimensional spacetime with compactly supported force, this would seem 
to definitively rule it out as a physically realistic equation of motion. However, 
it is an attractive conjecture that in this situation, all special solutions must be 
bounded. If so, then the space of physically relevant solutions coincides with 
the space of special solutions, and it has the physically expected dimension 3. 
Similar remarks apply to the DD equation in spacetimes of other dimensions. 

More generally, it is attractive to speculate that under reasonable hypotheses 
on a force not compactly supported (e.g., an exponentially decreasing force), 
all special solutions might be bounded. If this could be proved, it might be 
considered as a solution in principle to the problem posed in comment 3 of 
the Summary and Assessment section: Is there a generalized initial condition 
which will eliminate runaways in this situation? That is, we could identify 
"physical" (e.g., non-runaway) solutions with special solutions. However, since 
closed form special solutions are typically hard to come by, the practical problem 
of identifying the special solution in some constructive way would remain. 

The only thing which prevents us from directly applying the insights of 
the Ryabov/Driver theory to the DDI equation l|14l) is that the function sinh 
does not satisfy the required Lipschitz condition. But it seems possible that 
this may be merely a technical matter, and it seems reasonable to hope that 
the Ryabov/Driver theory might be extended to apply to this equation. If so, 
Theorem 5's conclusion that the rapidity 8(t) converges to a constant as r — > oo 
illustrates how the manifold of special solutions models the long-term behavior 
of the system. If every special solution converges to a constant in the future 
(as the theorem establishes for eventually vanishing force), then every solution 
must similarly converge to a constant (as the theorem also establishes). 

However, even if the Ryabov/Driver theory could be extended to cover the 
DD equation Q, there would be another difficulty which could be more fun- 
damental. This difficulty is that even for delay equations (|44H to which the 
Ryabov/Driver theory does apply, it only applies for sufficiently small t\. 

4 See, e.g., (Chicone, 2003) for an explicit bound ri < S which guarantees applicability 
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Unfortunately, in our application, it is questionable to assume that t\ is 
arbitrarily small. Indeed, a critic of an earlier version of the present work raised 
the following interesting objection. 

In Rohrlich's formulation, the classical charged particle is treated as a sphere 
of nonzero radius and the delay T\ is twice that radius (the time for light to 
traverse the sphere). If the radius is too small, then the classical electrostatic 
energy of such a sphere will exceed the typical energy required for the quantum 
effect of pair production, moving the problem out of the domain of classical 
electrodynamics and into the domain of quantum electrodynamics. But the DD 
equation is proposed solely as a classical equation of motion. 

Put differently, the considerations leading to the DD equation are not claimed 
to apply to arbitrarily small delays. The critic mistakenly thought that the 
present work assumed Caldirola's value of t\ — 4g 2 /(3mc 3 ) w 1.2 x 10 _23 sec, 
which he believed would be far less than could be the radius of any classical 
charged particle. 

Therefore, any proposed application of the Ryabov/Driver theory to the 
DD equation within the framework of motivation via a nonzero particle radius 
should address the question of whether the delay (equivalently, the particle's 
radius) is small enough for an extended Ryabov/Driver theory to apply. 

8 Appendix 2 

This appendix gives the details of the equivalence of 
da 1 

m 1 — = f(T)+m 2 [u\r-n)-u a (r-T 1 )u a (Ty(T)} , (6) 

and 

du 1 

m— = f + [u\T-\)~u a (T-\)u a (f)u\T)] . (10) 

We'll show that these equations are equivalent if the barred quantities are de- 
fined by: 

u(s) :— t\u{t\s) for all real numbers s, (46) 
f(s) :— ri/(ris)/m2 for all real numbers s, and (47) 
m := Tni/(rri2Ti) . (48) 

Call the time unit with respect to which |JBJ is written the "old time unit". 
First we change to a new time unit which equals t% old time units. If the old 
metric tensor is denoted tj := diag(l, —1,-1,-1), this can be accomplished by 
introducing a new metric fj := t^ 2 ?]. This changes the old proper time r to a 
new proper time 

f = t± t , T = Tlf 

of the Ryabov/Driver results. Unfortunately for our application, the bound 8 is inversely 
proportional to the Lipschitz constant, which is formally infinite for sinh (so that 8 is formally 
zero). 
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A worldline r i— > z(t) € R 4 parametrized by old proper time r corresponds 
to a "new" worldline f <—> z(f) parametrized by new proper time f, with z(f) = 
z(t), i.e., for any real number s, 

z(s) = z(tis) . 

The "new" worldline consists of the same set of points as the old worldline, but 
the paramctrization is different. 

The four- velocity u{?) of the new worldline is related to the old four- velocity 
u(t) by: 

dz 

u(f) :— — — tiu(t), equivalently, u(t) = t^ 1 u(t) 
ut 

Thus u and u are related by the simple transformations 



us 



tiu(t\s), and u(s) = t x u(t x s), for all real numbers s. (49) 



This gives a way to translate any statement about u into a corresponding state- 
ment about u, and vice versa. Similarly, 

du 2 du ... du _ 2 du , . 

— = n — , equivalently, — = t x —. (50) 

Next, divide JBJ by m,2, and substitute the above relations (|49|l and (|50|) . 
obtaining 

- rr 2 ^ = ™+T 1 - 1 [^(f -l)-u«(f -lK(fK(f)] . (51) 

Now inspection shows that after multiplication of l|51|) by n , I|1U|) will result if 
we define u, f, and m as in equations (|T7|) . and above. 



9 Appendix 3: Comment on the analysis of (Mo- 
niz and Sharp, 1977) of a non-relativistic ver- 
sion of the DD equation 

(Rohrlich, 1997) states the following, concerning the DD equation and a non- 
relativistic simplification of it which will be discussed below: 

"Returning to the overview of classical charged particle dynamics, 
one can summarize the present situation as very satisfactory: for a 
charged sphere there now exist equations of motion both relativisti- 
cally [this refers to the DD equation ©] and nonrelativistically that 
make sense and that are free of the problems that have plagued the 
theory for most of this century; these equations have no unphysical 
solution, no runaways, and no preaccelerations." 
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This resolves into the following claims: 

Claim 1: The DD equation has no preaccelerative solutions; 

Claim 2: The DD equation has no solutions which are runaway (in the future); 

Claim 3: The DD equation has no "unphysical solution". 

The analysis of this paper shows that all of these claims are at best optimistic 
and at worst false. 

We showed that Claim 1 is false as stated. However, we also noted that it 
can be reformulated (by adjoining appropriate generalized initial conditions) so 
as to become true for the special case of an field applied for a finite time. It is 
unknown whether the claim can be repaired for arbitrary fields, as discussed in 
Section 3. 

We showed that Claim 2 is true for a particle moving in one space dimension 
under a force which eventually vanishes; it is unknown whether it is true for 
general three-dimensional motion, as discussed in Appendix 1. Since the proof 
we gave for one dimension was nonroutine and special to that dimension, we 
suspect that new ideas will be required for a three-dimensional proof. 

The truth of Claim 3 depends on what one means by "unphysical" solutions, 
but if one considers both preacceleration and postacceleration "unphysical", 
then we showed that this claim cannot be repaired. 

The only evidence for any of these claims offered by (Rohrlich, 1997) or 
(Rohrlich, 1999) is an analysis by (Moniz and Sharp, 1977) of a non-relativistic 
version of the DD equation. Since the non-relativistic version is an entirely 
different equation, even if the claims were true for the non-relativistic version, 
they would not imply corresponding claims for the DD equation. However, a 
careful reading of (Moniz and Sharp, 1977) reveals that their analysis does not 
even prove the claims for their nonrelativistic equation. 

The following passage from (Moniz and Sharp, 1977) has sometimes been 
interpreted as making Claims 1 and 2 for their nonrelativistic version of the DD 
equation: 

[From the last paragraph of Section II, p. 2856] "Summarizing, we 
have found that including the effects of radiation reaction on a 
charged spherical shell results neither in runaway behavior nor in 
preacceleration if the charge radius of the shell L > cr . . ." 

This quote is from the part of the paper which treats only classical charged 
particles (as opposed to the quantum-mechanical treatment of later sections) . 

It is not clear whether Moniz and Sharp intended to assert that no solution 
can be either runaway or preaccelerative (i.e., Claims 1 and 2), or that given 
ordinary initial conditions specifying the position and velocity at a given time, 
there exists a solution satisfying these conditions and which is neither runaway 
nor preaccelerative. The aim of the present appendix is to dispel the confusion 
over these various claims by analyzing precisely what (Moniz and Sharp, 1977) 
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does prove. We shall see that what they actually show is closer to the latter 
than the former. 

The condition L > ct corresponds in our notation to the condition that 
our delay parameter n must not be too small. More specifically, the condition 
L > ct translates into a condition 

n > 5 , (52) 

where 8 is a certain positive parameter whose exact value will not be important 
to us. Henceforth we assume this condition. (Incidentally, for other values of 
the parameter t%, Moniz and Sharp do find unphysical solutions, either future 
runaway or oscillatory.) 

Their nonrelativistic version of the DD equation is: 

^ = h(t)-6[u(t-n)-u(t)] , (53) 

where u(i) represents the particle's three-dimensional velocity at time t, h is a 
force-like term (a three-dimensional force divided by certain constants), and b 
is a constant. This is their equation (2.10) on p. 2853 written in notation closer 
to ours. It can be obtained from the space part of the DD equation JSJ or 
by neglecting terms of quadratic or higher order in the velocity. 

It is clearly much simpler than the DD equation. In particular, it is what 
one might call a "linear" delay-differential equation, and this is critical to their 
analysis via Fourier transforms. 

Moniz and Sharp's analysis of the possibility of solutions which are runaway 
in the future assumes that h eventually vanishes, so that for sufficiently large 
times, it may be dropped from the right side of (|53|l . obtaining 

^ = 6[u(t-n)-u(*)] . (54) 

This equation admits exponential solutions 

u(t):=u Q e at , (55) 

with uo a constant vector and a a complex constant satisfying a — b[e~ aT1 — 1]. 
Moniz and Sharp then show that for t\ satisfying (|52J) . the real part of a is 
negative, so that these exponential solutions decay to as t becomes large. 

That is all they prove. But this only proves that no solutions of the simple 
exponential form (|55|l can be runaway in the future, which is not the same thing 
as proving that no solution can be runaway in the future. There are many 
solutions to (|54|l which are not exponential — indeed, the space of all solutions 
is infinite dimensional (as noted in our Section 3 and Appendix 1), whereas the 
space of exponential solutions only has dimension 3. Thus Moniz and Sharp's 
analysis proves neither that no solution of their equation (|53|l is runaway (Claim 
2) nor that there exists a non-runaway solution satisfying an arbitrary initial 
condition u(<o) = Uq. 
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(Moniz and Sharp, 1977) does give a prescription for writing down a formal 
Fourier transform of a solution of l|53|) for each forcing function h. (They do not 
address the question of whether the formal expression produced, which contains 
singularities, actually is the distributional Fourier transform of a solution.) This 
prescription produces solutions which are not preaccelerative (in the precise 
sense defined in Section 3). Thus they show (to the standard of rigor common 
in physics journals like Physical Review, not to mathematical standards) that 
there exist non- preaccelerative solutions. 

However, their proof is incomplete in one important respect — they do not 
show that there exists a nonpreaccelerative solution satisfying an arbitrary ini- 
tial condition u(io) = Uo. The nonpreaccelerative solution u which they produce 
is unique, leaving no room to satisfy general (ordinary) initial conditions. 

Nevertheless, it seems possible that this proof could be completed. Even if 
so, this would not prove Claim l's assertion that no solutions are preaccelerative, 
but they may not have intended to assert Claim 1. 

Indeed, there is an simple, explicit counterexample to Claim 1 for the non- 
relativistic DD equation i|53|) . Consider the case in which the force vanishes 
identically, i.e., h = in (J53J. For this case, any solution for which the accelera- 
tion du/dt does not vanish identically is preaccelerative — the particle accelerates 
before the force is applied (because the force is never applied). So, the expo- 
nential solutions u{t) = u e at previously noted in lj53|l are preaccelerative when 
Uo ^ 0, in contradiction to Claim 1. 

10 Appendix 4: Referees' reports; Part 1: His- 
tory) 

This paper was first submitted to the Journal of Mathematical Physics. The first 
referee's report was superficial and demonstrably incorrect in some important 
respects, so I requested a second referee. The second referee recommended 
rejection on the sole grounds that the problem which both (Rohrlich, 1999) and 
my paper addresses (the problem of finding a sensible equation of motion for a 
classical charged particle) is not "relevant for physics" . 

The first referee only objected to the methods of the paper, not to its rel- 
evance. The second referee, who was presumably furnished the first referee's 
objections to the paper's methods, did not question the paper's methods. 

I don't agree with the second referee's value judgment, but I can understand 
and respect it. His report was thoughtfully written and showed understanding 
of the paper. It is the only report of six (see below) which I can respect. 

I had submitted the paper to J. Math. Phys. because of its mathematical 
content, some of which is at a higher level than is typically published by journals 
like Physical Review. After the rejection by J. Math. Phys., I decided to submit 
it to Phys. Rev. D (PRD), which had published (Rohrlich, 1999). I thought 
(naively, it turned out), that having recently published Rohrlich's paper on the 
same topic, they would find it inconsistent to maintain that the problem which 
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both papers address was physically irrelevant. 

The story to follow may not always reflect well on the standards of PRD, 
so before giving it, I'd like to express my appreciation for the efficiency of their 
online submission procedure. The mechanics of dealing with them is a pleasure, 
compared to most journals. 

They are efficient and fast. Nobody likes a rejection, but an immediate 
rejection is infinitely more courteous than a rejection after several years of 
unanswered or inadequately answered correspondence. With other journals, 
the latter happens more often than it should. (I should also remark that J. 
Math. Phys. was also efficient.) 

Initially, PRD sent the paper to an anonymous referee and to Professor 
Rohrlich as an identified referee, the latter because the paper commented on 
his. Professor Rohrlich thought that the paper's analysis was invalid because it 
replaces his original equation © with the equivalent equation (JJJ) . I wrote him 
privately spelling out the equivalence (essentially the same analysis as Appendix 
2, which was added later), but he maintained his objection. 

The anonymous referee submitted a short, superficial report giving the im- 
pression that the paper merely quoted "existing theorems", and objecting that 
its subject was "interesting for the foundations of physics but not for phe- 
nomenology". (I wonder what he'd say about string theory!) 

I could understand this referee's value judgment if PRD were an experimental 
journal, but it seems inconsistent with its publication of Rohrlich's paper in 1999 
and the continued publication of papers on the same topic since that time. I 
wondered if this referee had actually read the paper (there was no internal 
evidence of that in the review), and I wondered if he might be looking for 
excuses to avoid dealing with its mathematics. 

Professor Rohrlich's objection to replacement of the equation © as he orig- 
inally wrote it with the equivalent version J7J) seemed so far out in left field that 
I felt I couldn't let it pass. So, I wrote the Editor explaining why I thought that 
Rohrlich had made a serious mistake, and was continuing to make it. 

He sent the paper, including Appendix 2 spelling out the equivalence of 
iJBJ and |0|, to a second anonymous referee (the "Third Referee", including 
Rohrlich). Third Referee reported that he "completely agree [s] with Professor 
Rohrlich's critique that 'results based on equation Q are physically meaningless 

I was flabbergasted, given that the issue in dispute is so elementary, and had 
been explained so carefully in Appendix 2. The report of Third Referee gave no 
indication that he had even read Appendix 2. Surely, if there were some error 
in the short Appendix 2, a conscientious referee should have pointed it out. 

By this time, I was so disgusted that it was hard not to simply walk away 
from the matter, but I eventually decided that for completeness and closure, I 
should protest this latest incompetence. So, I wrote the Editor, D. Nordstrom, 
explaining that I believed that, however unlikely it might seem, Third Referee's 
report was also seriously in error. 

I suggested the following way of resolving the matter to the satisfaction of all: 
submit the narrow issue of the equivalence of JJJJ and 10 to an identified referee 
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acceptable both to PRD and me. (I suggested Barry Simon, whom I have never 
met, but whose competence I respect, based on his published work.) I promised 
that if such a referee ruled that these two equations were not equivalent, I would 
immediately withdraw the paper and never submit another to PRD. 

He didn't reply to that, but he did send the paper back to Third Referee, 
who reversed himself, writing: 

"I have no doubt about the mathematical equivalence of equs. (6) 
and (7). (And I have read and confirmed the calculations in Ap- 
pendix 2!)" 

But this referee maintained his recommendation of rejection on the basis that 
(a) "the topic of a classical equation of motion does not belong to the most 
urgent questions in theoretical physics . . ." , and (b) the paper does not address 
the issue of the size of the postacceleration predicted by the DD equation. 

11 Appendix 4, Part 2: All of the referees' re- 
ports, in full 

This section of the previous version, titled as above, contained all of the referees' 
reports, along with discussions of them. The verbatim reports have been deleted 
from this version at the request of the arXiv. 

I am not convinced that it was legally necessary to remove the referees' 
reports. The legal issues involved seem to me obscure and unsettled. I shall not 
attempt to summarize them here. 

From a commonsense standpoint, it's hard to imagine why inclusion of the 
verbatim reports should be less acceptable than paraphrasing them (which is 
undisputedly allowed). If it should turn out that this is indeed what the law 
requires, then the law should be changed. The only reason for including the 
reports verbatim was to avoid any possible distortion of the referees' meanings. 

Some of the six reports were already paraphrased in the preceding section. 
To paraphrase them all and rewrite the discussions accordingly would require 
more work than I am willing to do. 

I will furnish the previous version, including the verbatim reports, to anyone 
interested, by email or other means. For further information, visit www.math.umb.edu/ 

12 Appendix 5: Solving the DD equation back- 
wards in time 

The referees for Physical Review D (PRD) raised the following objections (and 
only these) to the present paper. 

1. The version of the DD equation analyzed by the paper, (JJJ, is not physi- 
cally equivalent to the one proposed by Rohrlich, iJBJ. 
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After many months and many pages of careful explanation, the only inde- 
pendent referee to affirm this objection finally withdrew it. So, one hopes 
that this issue may be regarded as settled. 

2. The subject of the correct equation of motion for classical charged particles 
is of insufficient physical interest. 

3. The paper is too mathematical. 

4. The paper does not estimate the size of the postaccelerations predicted by 
the DD equation. 

Objection 2 seems odd, given that PRD recently published Rohrlich's paper 
|13| on precisely this topic (equations of motion for classical charged particles), 
along with papers on the same topic by other authors. However, taking it at face 
value, it occurred to me that it certainly shouldn't apply to a "Comment" paper, 
e.g., "Comment on Classical self-force by F. Rohrlich". (If it did apply, that 
would amount to an editorial declaration that the subject of Classical self-force 
was of insufficient physical interest to merit publication, not long after PRD had 
published it!) Also, such a paper, (which is required to be short and therefore 
could not contain a full mathematical analysis) would obviate Objection 3. 

I wrote to Editor D. Nordstrom asking if PRD would consider such a "Com- 
ment" paper, or if perhaps the editors had decided not to publish any form of 
the present work. After he didn't reply, I rewrote the paper as a Comment 
paper and submitted it on July 28, 2004. Its submission was acknowledged 
immediately, but by the end of December, I had heard nothing more. I wrote 
on December 28 enquiring about the status of the paper, and in particular if it 
had been sent to any referee. On January 28, 2005, I received a cryptic reply 
apologizing that the paper had been "misfiled" . and stating that it had been 
"sent out for review". (This is being written on January 31.) 

The presentation of the Comment paper is based on the fact that the solu- 
tion manifold for the DD equation is infinite- dimensional, rather than of finite 
dimension as expected (cf. Appendix 1). Thus to base a physically reason- 
able theory on this equation, one needs some auxiliary condition to reduce the 
solution manifold to finite dimension. 

That the DD equation's solution manifold is of infinite dimension is implicit 
in the present work and strongly suggested by the fact that the space of gen- 
eralized initial conditions is obviously infinite-dimensional. However, no careful 
proof that the solution manifold is infinite-dimensional was included above be- 
cause it was only used for motivation in peripheral discussions. 

The "Comment" paper contains a broad mathematical outline, but no sub- 
stantial proofs, for which the reader is referred to the present work. After 
submitting it, I realized that the present work didn't give a real proof that the 
solution manifold of the DD equation is infinite-dimensional. This appendix 
remedies this by furnishing a proof. 

Although the Comment paper uses Rohrlich's original form © of the DD 
equation, for notational simplicity, we use here the normalized DD equation Q , 
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so that the time units are normalized to make the delay t± equal to 1: 
dn i 

m- = f(r) + K(T-l)-/(T-lK(rK(r)] , (7) 

Recall that a generalized initial condition is a specification of u on some interval 
To — 1 < r < To of length 1. (For a general delay n, the interval of specification 
would be an interval of length t\.) Section 3 noted that to every such gen- 
eralized initial condition corresponds a unique solution defined on [tq , oo) . In 
other words, the DD equation can be solved forward in time starting with any 
generalized initial condition. 

We are now going to show how to solve the DD equation backward in time to 
obtain a unique solution on (— oo, To) for each generalized initial condition. In 
Section 3, this was done in a simple, ad hoc way for the DDI equation, which was 
all that was needed for the treatment given. But the Comment paper doesn't 
mention the DDI equation, so now we need to do it for the full DD equation 

o 

Combining the forward solution with the backward solution gives a unique 
solution on (— oo, oo) corresponding to each generalized initial condition on (to — 
1,to). We define the solution manifold to be the set of all solutions defined on 
(— oo, oo). 

Then it becomes obvious that the solution manifold is parametrized by the 
set of all generalized initial conditions on any given interval (to — 1, To) of length 
1. The technical issue is whether the assignment of a solution to each generalized 
initial condition is one-to-one (hence bijective). This is obvious if "solutions" 
are required to be defined on (— oo, oo) because then each generalized initial 
condition is the restriction of a solution. (It's not obvious if "solutions" are 
only required to be defined on (to, oo), though it can be proved.) 

To see how to solve the DD equation backward in time starting from a 
generalized initial condition, attempt to solve in J2J) for u(t— 1) in terms of m(t), 
du/dr, and f(r). This is impossible because the bracketed quantity is only the 
projection of u(t — 1) on u(r), but the equation does uniquely determine this 
projection. Hence it uniquely determines u(t — 1) up to an additive term f3u(r), 
(3 a scalar: 

dv l 

u{t - 1) = Pu{t) + [m— - f (t)] = 0u(t) + q(r) , 
where for brevity we introduce the name 

q(r) := [m^ - f (r)] 

for the bracketed quantity. 

Routine calculation reveals that the scalar (3 is uniquely determined by the 
necessary condition that the four- velocity u(t — 1) must be a future-pointing 
unit vector. The condition that u(t — 1) have unit norm gives a quadratic 
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equation for /?, and the condition that u(t — 1) be future-pointing singles out a 
unique solution to this quadratic equation. The result of this calculation is: 



u(t - 1) = [-u a (T)q a (r) + y/1 + (u»(T)q a (T)) 2 - (q a (r)q a (r)) 2 HT) + <z(r). 

Thus given the four-force /, the values of u(t) on an interval [to,to + 1], 
uniquely determine the values of u(t — 1) on this interval, Equivalently stated, 
the values of u(t) on [r , r + 1] uniquely determine the values of u{t) on [r — 
1,t ], and, by iteration, the values of u(t) on (oo,t ] 
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